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Abstract 

Molecular motors are essential components for the biophysical functions 
of the cell. Current quantitative understanding of how multiple motors move 
along a single track is not complete, even though models and theories for a single 
motor mechanochemistry abound. Recently, M.J.I. Miiller et al. have devel- 
oped a tug-of-war model to describe the bidirectional movement of the cargo 
(PNAS(2008) 105(12) P4609-4614). They found that the tug-of-war model ex- 
hibits several qualitative different motility regimes, which depend on the precise 
value of single motor parameters, and they suggested the sensitivity can be used 
by a cell to regulate its cargo traffic. In the present paper, we will carry out 
a further detailed theoretical analysis of the tug-of-war model. All the stable, 
i.e., biophysically observable, steady states and their stability domains can be 
obtained. Depending on values of the several parameters, tug-of-war model 
exhibits either uni-, bi- or tristability. In large motor number case, the steady 
state movement of the cargo, which is transported by two molecular motor 
species, is determined by the initial numbers of the motors which bound to the 



*School of Mathematical Sciences, Fudan University, Shanghai 200433, China 

' Shanghai Key Laboratory for Contemporary Applied Mathematics, Fudan University 

+ Centre for Computational Systems Biology, Fudan University (E-Mail: xyz@fudan.edu.cn) 



1 



track. For small motor number case, the movement of cargo may jump from 
one of the stable steady state to another. 

PACS: 87.16.Nn, 87.16.A-, 82.39.-k, 05.40.Jc 

Keywords: Tug-of-war, molecular motors, intracellular transport 

1 Introduction 

Molecular motors, including biological motor proteins such as kinesin [fl [2j [3j H] , 
dynein [5l [6] , mysion (TJ [HJ [9] and FqFi-ATP synthase [10] , are mechanochemical force 
generators which convert chemical or biochemical energy in the form of chemical po- 
tential into mechanical work in thermal environment [11] . The mechanochemical pro- 
cess is accomplished by individual macromolecules, immersed in an aqueous solution 
with the chemical potential, moving along a linear track. Many biological motor pro- 
teins move processively. For example, myosin slides along an actin filament, kinesin 
and dynein along microtubule (MT). All of them are adenosine triphosphate (ATP)- 
driven "directional walking machines" ( [T2"| [T3] ) : Kinesin moves towards the plus end 
of the MT and dynein towards the minus end. In comparison with the macroscopic 
engines driven by Carnot cycles, molecular motors have a high energy efficiency at 
about 50%, while the energy efficiency of a car is about 15%-20% [5j [HI [15]. Fur- 
thermore, the velocities of molecular motors are also fast with mean velocity be at 
about several hundreds nanometers per second [IB] . However, the most significant 
difference between the molecular motors and the macroscopic engines is that the for- 
mer are moving in a thermal noise dominated environment [IT]. So the movement of 
the molecular motors should be described stochastically, rather than determinately. 
Being able to convert and harvest energy with high efficiency on a mesoscopic scale 
makes molecular motors an exciting area of scientific research with potentially great 
innovative applications for energy production. 
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Great progress has been made in recent years in modeling the movement of molec- 
ular motors, including the mean field methods (THJ [Tj3 [H] , the Langevin stochastic 
dynamic methods [20], [21] and discrete stochastic methods [22] [23] [2U (25] [26] . How- 
ever, the existing models for a single molecular motor are not sufficient in predicting 
the recent experimental results: It is found that bidirectional motion of the cargo, 
which is carried by motor proteins, exhibits different patterns in different stages of em- 
bryonic development ( [27] ) . Following these recent experimental results ([28] [29] 50]). 
Lipowsky and his coworkers have developed the tug-of-war model for describing the 
movement of the cargo carried by processive motors, such as kinesin and dynein 
( [31] [32] l33l [3H [35]). In their model, the experimentally known single motor proper- 
ties are taken into account, so it is consistent with almost all experimental observa- 
tions and can make quantitative predictions for bidirectional transport of the cargo. 
Since cargo movement carried by a single motor protein via an elastic tether has 
been extensively studied in the past [36] [37], the focus of tug-of-war model is not on 
the detailed movement of cargo carried by a single motor per se, rather it concerns 
with the competition and cooperation of multiple motors on a single track (see the 
schematic depiction in Fig. [T]). 

In the present paper, we will give a further comprehensive mathematical analysis of 
tug-of-war model. Through detailed analysis, we find that the steady state movement 
of cargo is determined by the initial numbers of the two motor species which bound 
to the track of movement. Biophysically, the steady state is the only state that can 
be observed experimentally. At the same time, Monte Carlo simulations indicate the 
transition time from the initial state to the steady state is very short (see Figs. [7] 
[8]). Theoretically, the movement of the cargo has at most three stable steady states. 
If there exists two or three stable steady states, then many parameters of plus and 
minus motors have at least one critical point. The movement of cargo would change 
from one stable steady state to another if one of the parameters jumps from one 
side of its critical point to another side. In the following, we firstly introduce the 
tug-of-war model, and then give the detailed discussion gradually. 
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Figure 1: Schematic depiction of tug-of-war model: A cargo with N + = 3 plus motors 
(Kinesin) and i\L = 2 motors (Dynein) is pulled by a fluctuating number of motors 
bound to the microtubule. 

2 The tug-of-war model 

The tug-of-war model is developed by Reinhard Lipowsky's study group ([3T | |32 | 
1531 EH E5]) t° study the bidirectional transport of the cargo, in which the cargo is 
attached with N + plus and iV_ minus motors. Particularly, if N + = or iV_ = 0, it 
recovers the usual model for cooperate transport by a single motor species ([33] [38J ) . 
In this model, each motor species is characterized by six parameters, which can be 
measured in single molecular experiments (see Tab. [T]): (i) stall force F$ (pN) (ii) 
detachment force Fd (pN) (hi) unbinding rate eo (s _1 ) (iv) binding rate 7r (s _1 ) (v) 
forward velocity vf (/xm/s) and (vi) superstall velocity amplitude vb (nm/s). The 
motors bind to or unbind from a MT in a stochastic fashion, so that the cargo is 
pulled by n + < N + plus and n_ < iV„ minus motors, where n + and n_ fluctuate with 
time (see Fig. [Q. 

In tug-of-war model, it is assumed that, at every time t, the state of cargo with 
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Parameter 


Symbol 


Kinesin 1 


Dynein 


Stall force 


F s 


6pN 


l.lpN 


Detachment force 


F d 


3pN 


0.75pN 


Unbinding rate 


(Q 


Is' 1 


0.27s" 1 


Binding rate 




5s- 1 


1.6s" 1 


Forward velocity 


V F 


1/im/s 


0.65/zm/s 


superstall velocity amplitude 


VB 


6nm/s 


72nm/s 



Table 1: Single- motor parameters for kinesin 1 and cytoplasmic dynein ([31] and 
references therein). 

N + plus and N_ minus motors firmly attached to it is fully characterized by numbers 
n + and n_ of plus and minus motors that are bound to the MT. The state of cargo 
changes when a plus or a minus motor binds or unbinds to/from the MT (see Fig. 
[T]). The probability p(n + , n_, t) to have n + plus and n_ minus bound motors at time 
t can be described by the following Master equation: 

dp(n+,n-,t) r 

j t =\ N + - ("+ - IJF+pK- - l,n_,t) 

+ (n+ + l)e + (n + + 1, n_)p(n + + 1, n_, t) 
+ [N^ - (n_ - l)]7r_p(ra + ,n_ - 

+ (n_ + l)e+(n+, n_ + l)p(n + , + 1, t) (1) 

- [(JV+ - n + )n + + n + e + (n+, n_) 

+ (AL — r7,_)7r_ + n_e_(ra + , n_)]p{n +1 n_,t) 

1 <n + < N+ - 1 and 1 < n_ < iV_ - 1 

where 7r + (7r„) is the binding rate of a single plus (minus) motor to the MT, which 
depends only weakly on the load ([33]) and therefore is taken equal to zero-load 
binding rate 7ro+(7r -). e + (e_) is the unbinding rate of a single plus (minus) motor 
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from the MT, which increases exponentially with the applied force F: 

e±(F) = e ±exp(\F\/F d± ) (2) 

as measured for kinesin [39], where F^ is the detachment force. The governing 
equations for n + = 0,N + or n_ = 0,N_ are similar as (JTJ) except n + (N + ,n_) = 
-K-{n + , NJ) = and e+(0, n_) = e_(n+? 0) = 0. 

Under the assumptions that the motors act independently and feel each other only 
due to two effects: (i) opposing motors act as load, and (ii) identical motors share 
this load, Lipowsky and coworkers gave the following relation (see {34| ) 

n + F + = = F c (3) 

where F + (— F_) is the load felt by each plus (minus) motor. Eqs. d2J) ([3D imply 

e±(n+, n-) = e ± exp[|F c |/n±F d± ] (4) 

Here, the cargo force F c is determined by the condition that plus motors, which 
experience the force F c /n + , and minus motors, which experience the force — F c /n_, 
move with the same velocity v c , which is the cargo velocity: 

v c (n + ,n_) = v + (F c /n + ) = -v-(-F e /n+) (5) 

The same as in [31], the following piecewise linear force- velocity relation of a single 
motor is used in this paper: 

f v F (l - F/F s ) for F < F s 
v{F) = I (6) 

(v B (l-F/F s ) for F>F S 
where Vb is the absolute value of the superstall velocity amplitude, Vp is the zero- load 
forward velocity, F s is the stall force. 

3 The velocity of cargo and unbinding rates of mo- 
tors 

For the convenience of analysis in the following sections, we give the formulations 
of velocity of cargo and unbinding rates of plus and minus motors in this section. 
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(I) In case of "stronger plus motors", i.e. n + F s+ > n_F s _, Eqs. ((Sj) (JHD lead to the 
cargo force and velocity: 

VF+ + v B - 



FJn + ,n_, 

v F+ /n+F s+ + v B -/n^F s 

, \ n + F s+ - n_F s _ 



(7) 



n + F s+ /v F+ + n-F s -/v B ~ 

Using Eqs. PJ (JTj), the unbinding rates of plus and minus motors are: 

n T F s+ F s _(v F+ + v B _) 
[n + F a+ v B - + n„F s _v F+ )F d ± 



e±{n + ,n_) =e ±exp 
= : e ± exp 



,on+ + bnJ)F, 



d± 



where 



vb- , v F+ 



F s __ (v F+ + v B J) F s+ (v F+ + v B _] 

Let x = n + /N + , y = n_/N_ and c = N + /N_, then 



(9) 



e + (x,y) =e 0+ exp 
e-(x,y) =e _exp 



[acx + by)F M 
cx 



(10) 



(acx + by)F d . 

(II) In case of "stronger minus motors", i.e. n + F s+ < n_.F s _, the cargo force and 
velocity are: 

F c (n + ,n_) ■■ 



v B+ /n + F s+ + v F „/n_F s 



( \ n„F s ^-n + F s+ 

v c {n +J n_) = ^— ■ —- (11) 

n+F s+ v B+ + n-F s -/v F - 



yF s - — xcF t 



s+ 



xcF s+ /v B+ + yF s _/v F _ 
Similar as in (I), the unbinding rates of plus and minus motors are 



e+(x,y) =eo+exp 
e-(x,y) =e _exp 

in which 



y 



[acx + by)F d ^ 
cx 

[acx + by)F d . 



(12) 



VF- t V m 



F s _{v B+ + v F ^) F s+ (v B+ + v F J) 



(13) 
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The splitting boundary of case (I) and case (II) is n + F s+ = i.e. y = 

xcF s+ /F s _. 

(Ill) If an external force F ext is present, here F ext is taken to be positive if it points 
into the minus direction, then the force balance ([3]) becomes 



n + F + = -n-F+F, 



ext. 



In case of n + F s+ — F ext > n^F s _, carrying through the same calculation as for the 
case without external force leads to the velocity of cargo 

^K,n_)= n fy n f'-: F ;" (14) 
n + F s+ /v F+ + ri-F s ^/v B - 

The corresponding unbinding rates of the plus and minus motors are 

y + aF ext /N_ 



e + (x,y) =e +exp 
e-(x,y) =e _ exp 



(acx + by)F d+ 
cx - bF ext /N_ 



(15) 



(acx + by)F d _ 

(IV) If an external force F ext is present and n + F s+ — F ext < n_F s _, then the velocity 
of cargo is 

, . n + F s+ — n_F s _ — F ext 
v c (n+,n-) = — -— ■ -— (16) 

n + F s+ /v B + + n_F s ^/v F - 

and the unbinding rates of plus and minus motors are 

y + aF ext /N_ 



(acx + by)F ( 



d+ 



cx-bF ext /N^ U,) 



e+(x,y) =e +exp 
e-(x,y) =e -exp 

(acx + by)F d _ 

Similarly, the splitting boundary of case (III) and case (IV) is n + F s+ = n_F s _ + 
F ext , i.e. y = xcF s+ /F s „ - F ext /N-F s _. 

(V) More generally, if there exists an external force F ext and the friction coefficient 

of cargo is 7, then in the case of n + F s+ — F ext > n_F s _, the velocity of the cargo is 

n+F s + — n-F a - — F ext , s 

V + ; n+F a+ /v F+ + n-F a -/v B -+'y V ; 

and the unbinding rates of plus and minus motors are 

( , (y + a(F ext + 1Vc )/N\ 
e+(x,y) =e +exp , 

V (acx + by)F d+ J 

, s (cx-b(F ext + 1Vc )/N^ 
e-(x,y) =e _exp — — 

V (acx + by)F d 



On the other hand, if n + F s+ — F ext < n_F s _, then the velocity of cargo is 

/ , n + F s+ - n_F s _ - F ext 

v c (n +J n_) = — — 

n + F s+ /v B+ + n^F s _/v F ^ + 7 

and the unbinding rates of plus and minus motors are 

'y + a(F ext + jv c )/N_ 



e+(x,y) =e 0+ exp 
e-{x,y) =e _exp 



(acx + by)F d+ 
cx - b(F ext + jv c )/N. 



(acx + by)F d _ 

The splitting boundary of these two cases is also n + F s 
y = xcF s+ /F s _ - F ext /N„F S „. 



(20) 



(21) 



n-F s _ + F ext , i.e. 



4 The dynamics of motor numbers n + and n. 



For the sake of convenience, let 



r + = r + {n + , n_) :- 



s + = s + (n + , 7i- 



r_ = r_(n+, n_ 



(N+ - n + )u + 
n+e+(n +} n-) 
(N- - n_)7r_ 



(22) 



s_ = S-[n +1 n-) := n_e_^n + , n_J 

and A = r + + r_ + s + + s_. During time interval (t,t + At), the increase of plus motor 
number is 

•At 

A A 

In the limit At -> 0, ([23]) leads to 

dn. 



n + (t + At)-n + (t)=C-±-?±) I ' \ e ~ x dt= r -^ x S - 



(1 



— AAt 



) (23) 



r+ - s+ = (JV+ - n + )n+ - n+e + (n + , n_) 



(24) 



Similarly, the dynamics of minus motor number is 
(in. 



dt 



r + — s + = (N_ — n_)n_ — n_e_(n + , n_) 



(25) 
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Figure 2: The figures of functions f(x,y) = 0, g(x,y) = 0. The "+" ("-") means the 
function / (or g) is positive (negative) in the corresponding subdomains. 



So x = n + /N + , y = n_/iV_ satisfy 

dx 



- 7T+ - x[tt + + e+(x, y)\ := f(x, y) 

JL =7r _- y[ n _ + e _( X) y)\ ■= g( X} y) 



(26) 



As we all know, the steady state solutions (x*,y*) of the system (l26|) . which 
satisfy f(x*,y*) = and g(x*,y*) = 0, are stable if and only if the real parts of the 
two eigenvalues of the following matrix 



%(x*,y*) %(x*,y*) 



(27) 



are nonpositive. It is to say that 



dx 
df, 



, . ^dg , „ 



^-(x*,y*)<0 
dy 



df, , „dg, * „ N n 
^(^)^*)>0 



(2? 



To better understanding, the figures of functions f(x,y) = 0, g(x,y) = are 
plotted in figure [2J In view of conditions ( l28j) . to initial values xo = n + /N + , y = 
n_/iV_, if the point -Po^o, Ho) nes i n the subdomain I (II or III), then the final state is 
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Figure 3: The steady states of system (1261) . Where the unstable steady states are 
denoted by "*", the stable steady states are denoted by "®". If the initial state 
Po(xq, Do) lies in the subdomain I (II or III), then the final state is the stable steady 
state M i (M u or M 10 ).) 

stable steady state M i (M n or M 10 ) (see Fig. [3]). Theoretically, Um w 7^ 0, xm 01 ^ 0, 
but they are small than the accuracy of the numerical calculation used in this paper, 
so we simply regard them as 0. 

To further understand the properties of the stable steady state points, the figures 
of f(x, y) = and g(x, y) = with different values of parameters F s+ , Fd+, Fd-, 
vb+, vb-, Vp+, vf~, 7r + ,7r_, eo+,eo- and c = iV + /iV_ are plotted in Fig. H] and El El 
From the figures, one can find that system (f26|) might have one, two or three stable 
steady states, which depends on the values of the parameters. Given the initial value 
(xo,yo), the final steady state can be determined using the similar method as in Fig. 
|3] (Right). One can be easily know that, almost all of the parameters used in the 
tug-of-war model have one or two critical points, the final stable steady state would 
change if one of the parameters jumps from one side of its critical points to another 
side. 

Obviously, for iV + = or iV_ = (i.e. c = or c = oo), the tug-of-war model is 
reduced to the usual model for cooperate transport by a single motor species (minus 
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Figure 4: Figures of f(x,y) = 0,g(x,y) = for symmetric tug-of-war model, in 
which plus and minus motors have the same parameters, The unstable steady states 
are denoted by "*" , the stable steady states are denoted by "®" . 
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Figure 5: Asymmetric tug-of-war model: In this case, the system fl26l) might have 
one, two or three stable steady states. 




Figure 6: Tug-of-war model with external force F ext : In this case, the system (1261) 
might have two or three stable steady states. 
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Figure 7: For large motor numbers iV+, iV_ cases, the steady states is determined by 
the theoretical steady state n s + ps N + x s , nt ~ N^y s . Left: (n + (0)/N + , n_(0)/iV_) lie 
in subdomain (II); Middle: (n + (0)/N + ,n_(0)/N_) lie in subdomain (III); Right: 
(n + (0)/iV + ,n_(0)/iV_) lie in subdomain (I). 

or plus), and the only stable steady state is n + / (n + + eo+) for plus motor species or 
7r_/ (7r_ + e _) for minus motor species. The average velocity of the cargo at steady 
state is v c = v c (x*, 0) = vp+ if c = oo, and v c = v c (0, y*) = —vb- if c = 0, which are 
the velocities of a single motor. 

5 Comparison with Monte Carlo simulations 

Due to the above discussion, in large motor numbers limit iV + ,AL — > oo, the 
movement of the cargo might have one, two or three stable steady states. The final 
steady state is determined by the initial motor numbers n + (0) = N + x and n_(0) = 
N_y (see Fig. [7]). For example, in case of Fig. [3] (right), if (n+(0)/iV + , n_(0)/iV_) 
lies in subdomains (II) , the final steady state would be n s + m N + XMu,n s _ ~ ALy^n- 

However, if the numbers iV + , of molecular motors, which attached to the cargo, 
is finite or even small, the steady states numbers n s + and n s _ might be different with 
the theoretical values N + x* and N^y*. Theoretically, if Mi(xi,yi)(i = 1,2 or 3) are 
the stable steady points of the system (126]) . which can be regarded as the large motor 
numbers limit of fl24"l [25]) . then steady state numbers n s + and ni would lie in the 
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Figure 8: For small motor numbers N + ,N_, the final motor numbers n + ,n_ can 
change from one stable steady state to another. Left: The final motor numbers 
n + ,n^ change from N+xmu, ^-2/M n to N + xm w , N-Um 10 'i Right: The final motor 
numbers n +) n_ change from N + Xm 01 , N_y Moi to N+Xmw N_y Mll . 



neighborhoods of the theoretical values N+Xi and iV_yj. But, in small N + , iV_ cases, 
the steady state motor numbers n s + and n s _ can jump easily from the neighborhood 
of one of the theoretical stable steady state point (N + Xi, N-Ui) to the neighborhood 
of another theoretical stable steady state point (N + Xj, N_yj) (see Fig. [8]). For finite 
motor numbers N + , N_, the stepsize of the system fl26l) are Ax = 1/N+, Ay = 1/N_. 
So the smaller of motor numbers N + , AL, the easier for motor numbers n + , n_ to jump 
from one of the steady subdomains I, II or III to another. Intuitively, the probability 
that (n + /N + , n_/N_) lies in the neighborhood of the stable steady state point M, is 
proportional to the area of Mj's steady subdomain. Mathematically, the probability 
of motor numbers n+, n_ change from n^t , n_ to n+ , along trajectory S is 



n 

{Si,Si-\.\)&Sn 

n 

(S fc ,S fe+ i)GSj7 



n 



t 3 + 



(Sj,Sj+i)£SL 



n j+ + e j+ + TTj_ + Ci . 



7r fc+ + e fc+ + 7r fc _ + e fe . 



n 



(Si,Si + i)eSn 



(29) 
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where S L U S R U S v U S D = S, (Pi,P 2 ) G S R if and only if n+(P 2 ) = n+(Pi) + 
l,n_(P 2 ) = n_(Pi), (Pi,P 2 ) G S L if and only if n + (P 2 ) = n+(Pi) - l,n_(P 2 ) = 
n_(P t ), (P l5 P 2 ) G ^ if and only if n+(P 2 ) = n+(Pi),n_(P 2 ) = n_(P x ) + l, (P 1; P 2 ) G 
S'd if and only if n + (P 2 ) = n+(Pi), n_(P 2 ) = n_(Pi) — 1. So, theoretically, we can ob- 
tain the probability that motor numbers n+, n_ change from the neighborhood of one 
stable steady states to the neighborhood of another stable steady states. From these 
transition probabilities, we can know more details about the steady state movement 
of the cargo in this small iV + , iV_ cases. 

6 Conclusion and remarks 

In this paper, the steady state properties of the recent tug-of-war model, which 
is provided by Lipowsky et al to model the movement of cargo, which is transported 
by two motor species in the cell, is discussed. Biophysically, the stable steady states 
are the most important states, because the transition time to the stable steady state, 
as illustrated in this paper, is very short (see Figs. [7] and so almost all of the 
data are measured in stable steady states. Through the discussion in this paper, we 
can know that the final steady states of the movement of the cargo is determined by 
initial numbers of the plus and minus motors which are bounded to the microtubule. 
Certainly, the velocity and direction of the movement are also determined by other 
several parameters, such as N±, F s ±, n±, eo±, Fd±, vf±, vb±, Fext, 7- One can also find 
that, almost each of the parameters has critical points, which determine the stable 
steady velocity and direction of the cargo. It is most probable that, many of the 
parameters, including the numbers iV + and iV_ of plus and minus motors which are 
tightly attached to the cargo, and the initial binding numbers n + (0) and n_(0), can 
be determined by the biochemical environment and properties of the cargos, so some 
of which can be transported from the plus end to the minus end, and others can be 
transported reversely. 

Acknowledgments 



16 



This work was funded by National Natural Science Foundation of China (Grant 
No. 10701029). The author thanks professor Hong Qian of University of Washington 
for his help to complete this research. The author also thanks the reviewers for their 
help to improve the quality of this paper. 

References 

[1] David D. Hackney. Processive motor movement. Science, 316:58-59, 2007. 

[2] N. J. Carter and R. A. Cross. Mechanics of the kinesin step. Nature, 435:308-312, 
2005. 

[3] Y. Taniguchi, M. Nishiyama, Y. Ishhi, and T. Yanagida. Entropy rectifies the 
brownian step of kinesin. Nature Chemical Biology, 1:342-347, 2005. 

[4] Yunxin Zhang. Three phase model of the processive motor protein kinesin. Bio- 
physical Chemistry, 136:19-22, 2008. 

[5] R. D. Vale. The molecular motor toolbox for intracellular transport. Cell, 
112:467-480, 2003. 

[6] H. Sakakibara, H. Kojima, Y. Sakai, E. Katayama, and K. Oiwa. Inner-arm 
dynein c of chlamydomonas flagella is a single-headed processive motor. Nature, 
400:596-589, 1999. 

[7] A. M. Hooft, E. J. Maki, K. K. Cox, and J. E. Baker. An accelerated state of 
myosin-based actin motility. Biochemistry, 46:3513-3520, 2007. 

[8] J. Christof, M. Gebhardt, Anabel E.-M. Clemen, Johann Jaud, and Matthias 
Rief. Myosin- v is a mechanical ratchet. Proceedings of the National Academy of 
Sciences of the United States of America, 103:8680-8685, 2006. 

[9] Katsuyuki Shiroguchi and Kazuhiko Kinosita Jr. Myosin v walks by lever brow- 
nian motion. Science, 316:1208-1212, 2007. 



17 



[10] H. Noji, R. Yasuda, M. Yoshida, and Jr. K. Kinosita. Direct observation of the 
rotation of fl-atpase. Nature, 386:299-302, 1997. 

[11] J. Howard. Mechanics of Motor Proteins and the Cytoskeleton. Sinauer Asso- 
ciates, Sunderland, MA, 2001. 

[12] Hongyun Wang and George Oster. Energy transduction in the fl motor ofatp 
synthase. Nature, 396:279-280, 1998. 

[13] Masayoshi Nishiyama, Hideo Higuchi, and Toshio Yanagida. Chemo mechanical 
coupling of the forward and backward steps of single kinesin molecules. Nature 
Cell Biology, 4:790-797, 2002. 

[14] Hongyun Wang. Chemical and mechanical efficiencies of molecular motors and 
implications for motor mechanisms. Journal of Physics: Condensed Matter, 
17:S3997-S4014, 2005. 

[15] Yunxin Zhang. The efficiency of the molecular motors. Journal of Statistical 
Physics, 134:669-679, 2009. 

[16] K. Svoboda and S.M. Block. Force and velocity measured for single kinesin 
molecules. Cell, 77:773-784, 1994. 

[17] Peter Reimann. Brownian motors: noisy transport far from equilibrium. Physics 
Reports, 361:57, 2002. 

[18] Hongyun Wang. Several issues in modeling molecular motors. Journal of Com- 
putational and Theoretical Nanoscience, 5:1-35, 2008. 

[19] Hong Qian. The mathematical theory of molecular motor movement and chemo- 
mechanical energy transduction. Journal of Mathematical Chemistry, 27:219- 
234, 2000. 

[20] P. Reimann, C. Van den Broeck, P. Hanggi H. Linke, J.M. Rubi, and A. Perez- 
Madrid. Giant acceleration of free diffusion by use of tilted periodic potentials. 
Physical Review Letters, 87:010602, 2001. 



[21] R.D. Astumian. Thermodynamics and kinetics of a brownian motor. Science, 
276:917-922, 1997. 

[22] R Dean Astumian. Biasing the random walk of a molecular motor. Journal of 
Physics: Condensed Matter, 17:S3753-S3766, 2005. 

[23] Steffen Liepelt and Reinhard Lipowsky. Kinesins network of chemomechanical 
motor cycles. Physical Review Letters, 98(25) :258102, 2007. 

[24] Michael E. Fisher and Anatoly B. Kolomeisky. Simple mechanochemistry de- 
scribes the dynamics of kinesin molecules. Proceedings of the National Academy 
of Sciences of the United States of America, 98(14):7748-7753, 2001. 

[25] B. Derrida, J. L. Lebowitz, and E. R. Speer. Shock profiles for the asymmetric 
simple exclusion process in one dimension. Journal of Statistical Physics, 89:135- 
167, 1997. 

[26] Hong Qian. A simple theory of motor protein kinetics and energetics. Biophysical 
Chemistry, 67:263-267, 1997. 

[27] Steven P. Gross, Michael A. Welte, Steven M. Block, and Eric F. Wieschaus. Co- 
ordination of opposite-polarity microtubule motors. The Journal of Cell Biology, 
156:715-724, 2002. 

[28] Sean W. Deacon, Anna S. Serpinskaya, Patricia S. Vaughan, Monica Lopez Fa- 
narraga, Isabelle Vernos, Kevin T. Vaughan, and Vladimir I. Gelfand. Dynactin 
is required for bidirectional organelle transport. The Journal of Cell Biology, 
160:297-301, 2003. 

[29] Michael A. Welte, Silvia Cermelli, John Griner, Arturo Viera, Yi Guo, Dae- 
Hwan Kim, Joseph G. Gindhart, and Steven P. Gross. Regulation of lipid-droplet 
transport by the perilipin homolog lsd2. Current Biology, 15:1266-1275, 2005. 

[30] G. A. Smith, L. Pomeranz S. P. Gross, and L. W. Enquist. Local modulation of 
plus-end transport targets herpesvirus entry and egress in sensory axons. Pro- 

19 



ceedings of the National Academy of Sciences of the United States of America, 
101(45):16034-16039, 2004. 

[31] Melanie J.I. Muller, Stefan Klumpp, and Reinhard Lipowsky. Tug-of-war as 
a cooperative mechanism for bidirectional cargo transport by molecular motors. 
Proceedings of the National Academy of Sciences of the United States of America, 
105(12):4609-4614, 2008. 

[32] Janina Beeg, Stefan Klumpp, Rumiana Dimova, Ruben Serral Gracia, Eberhard 
Unger, , and Reinhard Lipowsky. Transport of beads by several kinesin motors. 
Biophysical Journal, 94:532-541, 2008. 

[33] Stefan Klumpp and Reinhard Lipowsky. Cooperative cargo transport by several 
molecular motors. Proceedings of the National Academy of Sciences of the United 
States of America, 102:17284-17289, 2005. 

[34] Melanie J.I. Muller, Janina Beeg, Rumiana Dimova, Stefan Klumpp, and Rein- 



hard Lipowsky. Traffic by small teams of molecular motors. arXiv:0807.0964vl 
[cond-mat.stat-mech]. 

[35] Melanie J.I. Muller, Stefan Klumpp, and Reinhard Lipowsky. Motility states 
of molecular motors engaged in a stochastic tug-of-war. Journal of Statistical 
Physics, 133:1059-1081, 2008. 

[36] Yi der Chen, Bo Yan, and Robert J. Rubin. Fluctuations and randomness of 
movement of the bead powered by a single kinesin molecule in a force-clamped 
motility assay: Monte carlo simulations. Biophysical Journal, 83:2360-2369, 
2002. 

[37] Timothy C. Elston and Charles S. Peskin. The role of protein flexibility in 
molecular motor function: Coupled diffusion in a tilted periodic potential. SIAM 
Journal on Applied Mathematics, 60(3):842-867, 2000. 

[38] Frank Julicher and Jacques Prost. Cooperative molecular motors. Physical Re- 
view Letters, 75:2618-2621, 1995. 

20 



K. Svoboda and S.M. Block. Force and velocity measured for single kinesin 
molecules. Cell, 77:773-784, 1994. 



21 



Properties of tug-of-war model for cargo 
transport by molecular motors 

Yunxin Zhang*T+ 



Abstract 

Molecular motors are essential components for the biophysical functions 
of the cell. Current quantitative understanding of how multiple motors move 
along a single track is not complete, even though models and theories for a single 
motor mechanochemistry abound. Recently, M.J.I. Miiller et al. have devel- 
oped a tug-of-war model to describe the bidirectional movement of the cargo 
(PNAS(2008) 105(12) P4609-4614). They found that the tug-of-war model ex- 
hibits several qualitative different motility regimes, which depend on the precise 
value of single motor parameters, and they suggested the sensitivity can be used 
by a cell to regulate its cargo traffic. In the present paper, we will carry out 
a further detailed theoretical analysis of the tug-of-war model. All the stable, 
i.e., biophysically observable, steady states and their stability domains can be 
obtained. Depending on values of the several parameters, tug-of-war model 
exhibits either uni-, bi- or tristability. In large motor number case, the steady 
state movement of the cargo, which is transported by two molecular motor 
species, is determined by the initial numbers of the motors which bound to the 
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track. For small motor number case, the movement of cargo may jump from 
one of the stable steady state to another. 

PACS: 87.16.Nn, 87.16.A-, 82.39.-k, 05.40.Jc 

Keywords: Tug-of-war, molecular motors, intracellular transport 

1 Introduction 

Molecular motors, including biological motor proteins such as kinesin [H, [2J, [31 ?] , 
dynein [U,[5], mysion [E1I3IH] and F Fi-ATP synthase [9j, are mechanochemical force 
generators which convert chemical or biochemical energy in the form of chemical 
potential into mechanical work in thermal environment [10] . The mechanochemical 
process is accomplished by individual macromolecules, immersed in an aqueous solu- 
tion with the chemical potential, moving along a linear track. Many biological motor 
proteins move processively. For example, myosin slides along an actin filament, ki- 
nesin and dynein along microtubule (MT). All of them are adenosine triphosphate 
(ATP)-driven "directional walking machines" ([IUEI2]): Kinesin moves towards the 
plus end of the MT and dynein towards the minus end. In comparison with the macro- 
scopic engines driven by Carnot cycles, molecular motors have a high energy efficiency 
at about 50%, while the energy efficiency of a car is about 15%-20% [H [131 EH]. Fur- 
thermore, the velocities of molecular motors are also fast with mean velocity be at 
about several hundreds nanometers per second [T5] . However, the most significant 
difference between the molecular motors and the macroscopic engines is that the for- 
mer are moving in a thermal noise dominated environment [16]. So the movement of 
the molecular motors should be described stochastically, rather than determinately. 
Being able to convert and harvest energy with high efficiency on a mesoscopic scale 
makes molecular motors an exciting area of scientific research with potentially great 
innovative applications for energy production. 
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Great progress has been made in recent years in modeling the movement of molec- 
ular motors, including the mean field methods [TTJ [TH1 OH] , the Langevin stochastic 
dynamic methods [191 [20] an d discrete stochastic methods [211 (22J [23J [2U ?]■ However, 
the existing models for a single molecular motor are not sufficient in predicting the 
recent experimental results: It is found that bidirectional motion of the cargo, which 
is carried by motor proteins, exhibits different patterns in different stages of embry- 
onic development ( [25] ) . Following these recent experimental results ([261 123 25]). 
Lipowsky and his coworkers have developed the tug-of-war model for describing the 
movement of the cargo carried by processive motors, such as kinesin and dynein 
( [29[ 1301 I3TI [321 ?]). In their model, the experimentally known single motor proper- 
ties are taken into account, so it is consistent with almost all experimental observa- 
tions and can make quantitative predictions for bidirectional transport of the cargo. 
Since cargo movement carried by a single motor protein via an elastic tether has 
been extensively studied in the past [?, ?], the focus of tug-of-war model is not on 
the detailed movement of cargo carried by a single motor per se, rather it concerns 
with the competition and cooperation of multiple motors on a single track (see the 
schematic depiction in Fig. [T]). 

In the present paper, we will give a further comprehensive mathematical analysis of 
tug-of-war model. Through detailed analysis, we find that the steady state movement 
of cargo is determined by the initial numbers of the two motor species which bound 
to the track of movement. Biophysically, the steady state is the only state that can 
be observed experimentally. At the same time, Monte Carlo simulations indicate the 
transition time from the initial state to the steady state is very short (see Figs. 
[HJ. Theoretically, the movement of the cargo has at most three stable steady states. 
If there exists two or three stable steady states, then many parameters of plus and 
minus motors have at least one critical point. The movement of cargo would change 
from one stable steady state to another if one of the parameters jumps from one 
side of its critical point to another side. In the following, we firstly introduce the 
tug-of-war model, and then give the detailed discussion gradually. 
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Figure 1: Schematic depiction of tug-of-war model: A cargo with N + = 3 plus motors 
(Kinesin) and i\L = 2 motors (Dynein) is pulled by a fluctuating number of motors 
bound to the microtubule. 

2 The tug-of-war model 

The tug-of-war model is developed by Reinhard Lipowsky's study group ([29 | [30 | 
EH E21 ?]) to study the bidirectional transport of the cargo, in which the cargo is 
attached with N + plus and iV_ minus motors. Particularly, if N + = or iV_ = 0, it 
recovers the usual model for cooperate transport by a single motor species ([31] |33j). 
In this model, each motor species is characterized by six parameters, which can be 
measured in single molecular experiments (see Tab. [T]): (i) stall force F$ (pN) (ii) 
detachment force Fa (pN) (iii) unbinding rate eo (s _1 ) (iv) binding rate ttq (s _1 ) (v) 
forward velocity vf (/xm/s) and (vi) superstall velocity amplitude vb (nm/s). The 
motors bind to or unbind from a MT in a stochastic fashion, so that the cargo is 
pulled by n + < N + plus and n_ < iV„ minus motors, where n + and n_ fluctuate with 
time (see Fig. [Q. 

In tug-of-war model, it is assumed that, at every time t, the state of cargo with 
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Parameter 


Symbol 


Kinesin 1 


Dynein 


Stall force 


F s 


6pN 


l.lpN 


Detachment force 


F d 


3pN 


0.75pN 


Unbinding rate 


(Q 


Is' 1 


0.27s" 1 


Binding rate 




5s- 1 


1.6s" 1 


Forward velocity 


V F 


1/im/s 


0.65/zm/s 


superstall velocity amplitude 


VB 


6nm/s 


72nm/s 



Table 1: Single- motor parameters for kinesin 1 and cytoplasmic dynein ([29] and 
references therein). 

N + plus and N_ minus motors firmly attached to it is fully characterized by numbers 
n + and n_ of plus and minus motors that are bound to the MT. The state of cargo 
changes when a plus or a minus motor binds or unbinds to/from the MT (see Fig. 
[T]). The probability p(n + , n_, t) to have n + plus and n_ minus bound motors at time 
t can be described by the following Master equation: 

dp(n+,n-,t) r 

j t =\ N + - ("+ - IJF+pK- - l,n_,t) 

+ (n+ + l)e + (n + + 1, n_)p(n + + 1, n_, t) 
+ [N^ - (n_ - l)]7r_p(ra + ,n_ - 

+ (n_ + l)e+(n+, n_ + n_ + 1, t) (1) 

- [(JV+ - n + )n + + n + e + (n+, n_) 

+ (AL — r7,_)7r_ + n_e_(ra + , n_)]p{n +1 n_,t) 

1 <n + < N+ - 1 and 1 < n_ < iV_ - 1 

where 7r + (7r„) is the binding rate of a single plus (minus) motor to the MT, which 
depends only weakly on the load ([SI]) and therefore is taken equal to zero-load 
binding rate 7ro+(7r -). e + (e_) is the unbinding rate of a single plus (minus) motor 
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from the MT, which increases exponentially with the applied force F: 

e±(F) = e ±exp(\F\/F d± ) (2) 

as measured for kinesin [31], where F^ is the detachment force. The governing 
equations for n + = 0,N + or n_ = 0,N_ are similar as (JTJ) except n + (N + ,n_) = 
-K-{n + , NJ) = and e+(0, n_) = e_(n+? 0) = 0. 

Under the assumptions that the motors act independently and feel each other only 
due to two effects: (i) opposing motors act as load, and (ii) identical motors share 
this load, Lipowsky and coworkers gave the following relation (see [32J) 

n + F + = = F c (3) 

where F + (— F_) is the load felt by each plus (minus) motor. Eqs. d2J) ([3D imply 

e±(n+, n-) = e ± exp[|F c |/n±F d± ] (4) 

Here, the cargo force F c is determined by the condition that plus motors, which 
experience the force F c /n + , and minus motors, which experience the force — F c /n_, 
move with the same velocity v c , which is the cargo velocity: 

v c (n + ,n_) = v + (F c /n + ) = -v-(-F e /n+) (5) 

The same as in [29J, the following piecewise linear force- velocity relation of a single 
motor is used in this paper: 

f v F (l - F/F s ) for F < F s 
v{F) = I (6) 

(v B (l-F/F s ) for F>F S 
where Vb is the absolute value of the superstall velocity amplitude, Vp is the zero- load 
forward velocity, F s is the stall force. 

3 The velocity of cargo and unbinding rates of mo- 
tors 

For the convenience of analysis in the following sections, we give the formulations 
of velocity of cargo and unbinding rates of plus and minus motors in this section. 
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(I) In case of "stronger plus motors", i.e. n + F s+ > n_F s _, Eqs. ((Sj) (JHD lead to the 
cargo force and velocity: 

VF+ + v B - 



FJn + ,n_, 

v F+ /n+F s+ + v B -/n^F s 

, \ n + F s+ - n_F s _ 



(7) 



n + F s+ /v F+ + n-F s -/v B ~ 

Using Eqs. PJ (JTj), the unbinding rates of plus and minus motors are: 

n T F s+ F s _(v F+ + v B _) 
[n + F a+ v B - + n„F s _v F+ )F d ± 



e±{n + ,n_) =e ±exp 
= : e ± exp 



,on+ + bnJ)F, 



d± 



where 



vb- , v F+ 



F s __ (v F+ + v B J) F s+ (v F+ + v B _] 

Let x = n + /N + , y = n_/N_ and c = N + /N_, then 



(9) 



e + (x,y) =e 0+ exp 
e-(x,y) =e _exp 



[acx + by)F M 
cx 



(10) 



(acx + by)F d . 

(II) In case of "stronger minus motors", i.e. n + F s+ < n_.F s _, the cargo force and 
velocity are: 

F c (n + ,n_) ■■ 



v B+ /n + F s+ + v F „/n_F s 



( \ n„F s ^-n + F s+ 

v c {n +J n_) = ^— ■ —- (11) 

n+F s+ v B+ + n-F s -/v F - 



yF s - — xcF t 



s+ 



xcF s+ /v B+ + yF s _/v F _ 
Similar as in (I), the unbinding rates of plus and minus motors are 



e+(x,y) =eo+exp 
e-(x,y) =e _exp 

in which 



y 



[acx + by)F d ^ 
cx 

[acx + by)F d . 



(12) 



VF- t V m 



F s _{v B+ + v F ^) F s+ (v B+ + v F J) 



(13) 
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The splitting boundary of case (I) and case (II) is n + F s+ = i.e. y = 

xcF s+ /F s _. 

(Ill) If an external force F ext is present, here F ext is taken to be positive if it points 
into the minus direction, then the force balance ([3]) becomes 



n + F + = -n-F+F, 



ext. 



In case of n + F s+ — F ext > n^F s _, carrying through the same calculation as for the 
case without external force leads to the velocity of cargo 

^K,n_)= n fy n f'-: F ;" (14) 
n + F s+ /v F+ + ri-F s ^/v B - 

The corresponding unbinding rates of the plus and minus motors are 

y + aF ext /N_ 



e + (x,y) =e +exp 
e-(x,y) =e _ exp 



(acx + by)F d+ 
cx - bF ext /N_ 



(15) 



(acx + by)F d _ 

(IV) If an external force F ext is present and n + F s+ — F ext < n_F s _, then the velocity 
of cargo is 

, . n + F s+ — n_F s _ — F ext 
v c (n+,n-) = — -— ■ -— (16) 

n + F s+ /v B + + n_F s ^/v F - 

and the unbinding rates of plus and minus motors are 

y + aF ext /N_ 



(acx + by)F ( 



d+ 



cx-bF ext /N^ U,) 



e+(x,y) =e +exp 
e-(x,y) =e -exp 

(acx + by)F d _ 

Similarly, the splitting boundary of case (III) and case (IV) is n + F s+ = n_F s _ + 
F ext , i.e. y = xcF s+ /F s „ - F ext /N-F s _. 

(V) More generally, if there exists an external force F ext and the friction coefficient 

of cargo is 7, then in the case of n + F s+ — F ext > n_F s _, the velocity of the cargo is 

n+F s + — n-F a - — F ext , s 

V + ; n+F a+ /v F+ + n-F a -/v B -+'y V ; 

and the unbinding rates of plus and minus motors are 

( , (y + a(F ext + 1Vc )/N\ 
e+(x,y) =e +exp , 

V (acx + by)F d+ J 

, s (cx-b(F ext + 1Vc )/N^ 
e-(x,y) =e _exp — — 

V (acx + by)F d 



On the other hand, if n + F s+ — F ext < n_F s _, then the velocity of cargo is 

/ , n + F s+ - n_F s _ - F ext 

v c (n +J n_) = — — 

n + F s+ /v B+ + n^F s _/v F ^ + 7 

and the unbinding rates of plus and minus motors are 

'y + a(F ext + jv c )/N_ 



e+(x,y) =e 0+ exp 
e-{x,y) =e _exp 



(acx + by)F d+ 
cx - b(F ext + jv c )/N. 



(acx + by)F d _ 

The splitting boundary of these two cases is also n + F s 
y = xcF s+ /F s _ - F ext /N„F S „. 



(20) 



(21) 



n-F s _ + F ext , i.e. 



4 The dynamics of motor numbers n + and n. 



For the sake of convenience, let 



r + = r + {n + , n_) :- 



s + = s + (n + , 7i- 



r_ = r_(n+, n_ 



(N+ - n + )u + 
n+e+(n +} n-) 
(N- - n_)7r_ 



(22) 



s_ = S-[n +1 n-) := n_e_^n + , n_J 

and A = r + + r_ + s + + s_. During time interval (t,t + At), the increase of plus motor 
number is 

•At 

A A 

In the limit At -> 0, ([23]) leads to 

dn. 



n + (t + At)-n + (t)=C-±-?±) I ' \ e ~ x dt= r -^ x S - 



(1 



— AAt 



) (23) 



r+ - s+ = (JV+ - n + )n+ - n+e + (n + , n_) 



(24) 



Similarly, the dynamics of minus motor number is 
(in. 



dt 



r + — s + = (N_ — n_)n_ — n_e_(n + , n_) 



(25) 
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Figure 2: The figures of functions f(x,y) = 0, g(x,y) = 0. The "+" ("-") means the 
function / (or g) is positive (negative) in the corresponding subdomains. 



So x = n + /N + , y = n_/iV_ satisfy 

dx 



- 7T+ - x[tt + + e+(x, y)\ := f(x, y) 

JL =7r _- y[ n _ + e _( X) y)\ ■= g( X} y) 



(26) 



As we all know, the steady state solutions (x*,y*) of the system (l26|) . which 
satisfy f(x*,y*) = and g(x*,y*) = 0, are stable if and only if the real parts of the 
two eigenvalues of the following matrix 



%(x*,y*) %(x*,y*) 



(27) 



are nonpositive. It is to say that 



dx 
df, 



, . ^dg , „ 



^-(x*,y*)<0 
dy 



df, , „dg, * „ N n 
^(^)^*)>0 



(2? 



To better understanding, the figures of functions f(x,y) = 0, g(x,y) = are 
plotted in figure [2J In view of conditions ( l28j) . to initial values xo = n + /N + , y = 
n_/iV_, if the point -Po^o, Ho) nes i n the subdomain I (II or III), then the final state is 
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Figure 3: The steady states of system (1261) . Where the unstable steady states are 
denoted by "*", the stable steady states are denoted by "®". If the initial state 
Po(xq, Do) lies in the subdomain I (II or III), then the final state is the stable steady 
state M i (M u or M 10 ).) 

stable steady state M i (M n or M 10 ) (see Fig. [3]). Theoretically, Um w 7^ 0, xm 01 ^ 0, 
but they are small than the accuracy of the numerical calculation used in this paper, 
so we simply regard them as 0. 

To further understand the properties of the stable steady state points, the figures 
of f(x, y) = and g(x, y) = with different values of parameters F s+ , Fd+, Fd-, 
vb+, vb-, Vp+, vf~, 7r + ,7r_, eo+,eo- and c = iV + /iV_ are plotted in Fig. H] and El El 
From the figures, one can find that system (f26|) might have one, two or three stable 
steady states, which depends on the values of the parameters. Given the initial value 
(xo,yo), the final steady state can be determined using the similar method as in Fig. 
|3] (Right). One can be easily know that, almost all of the parameters used in the 
tug-of-war model have one or two critical points, the final stable steady state would 
change if one of the parameters jumps from one side of its critical points to another 
side. 

Obviously, for iV + = or iV_ = (i.e. c = or c = oo), the tug-of-war model is 
reduced to the usual model for cooperate transport by a single motor species (minus 
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Figure 4: Figures of f(x,y) = 0,g(x,y) = for symmetric tug-of-war model, in 
which plus and minus motors have the same parameters, The unstable steady states 
are denoted by "*" , the stable steady states are denoted by "®" . 



12 



Figure 5: Asymmetric tug-of-war model: In this case, the system fl26l) might have 
one, two or three stable steady states. 




Figure 6: Tug-of-war model with external force F ext : In this case, the system (1261) 
might have two or three stable steady states. 
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Figure 7: For large motor numbers iV+, iV_ cases, the steady states is determined by 
the theoretical steady state n s + ps N + x s , nt ~ N^y s . Left: (n + (0)/N + , n_(0)/iV_) lie 
in subdomain (II); Middle: (n + (0)/N + ,n_(0)/N_) lie in subdomain (III); Right: 
(n + (0)/iV + ,n_(0)/iV_) lie in subdomain (I). 

or plus), and the only stable steady state is n + / (n + + eo+) for plus motor species or 
7r_/ (7r_ + e _) for minus motor species. The average velocity of the cargo at steady 
state is v c = v c (x*, 0) = vp+ if c = oo, and v c = v c (0, y*) = —vb- if c = 0, which are 
the velocities of a single motor. 

5 Comparison with Monte Carlo simulations 

Due to the above discussion, in large motor numbers limit iV + ,AL — > oo, the 
movement of the cargo might have one, two or three stable steady states. The final 
steady state is determined by the initial motor numbers n + (0) = N + x and n_(0) = 
N_y (see Fig. [7]). For example, in case of Fig. [3] (right), if (n+(0)/iV + , n_(0)/iV_) 
lies in subdomains (II) , the final steady state would be n s + m N + XMu,n s _ ~ ALy^n- 

However, if the numbers iV + , of molecular motors, which attached to the cargo, 
is finite or even small, the steady states numbers n s + and n s _ might be different with 
the theoretical values N + x* and N^y*. Theoretically, if Mi(xi,yi)(i = 1,2 or 3) are 
the stable steady points of the system (126]) . which can be regarded as the large motor 
numbers limit of fl24"l [25]) . then steady state numbers n s + and ni would lie in the 
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Figure 8: For small motor numbers N + ,N_, the final motor numbers n + ,n_ can 
change from one stable steady state to another. Left: The final motor numbers 
n + ,n^ change from N+xmu, ^-2/M n to N + xm w , N-Um 10 'i Right: The final motor 
numbers n +) n_ change from N + Xm 01 , N_y Moi to N+Xmw N_y Mll . 



neighborhoods of the theoretical values N+Xi and iV_yj. But, in small N + , iV_ cases, 
the steady state motor numbers n s + and n s _ can jump easily from the neighborhood 
of one of the theoretical stable steady state point (N + Xi, N-Ui) to the neighborhood 
of another theoretical stable steady state point (N + Xj, N_yj) (see Fig. [8]). For finite 
motor numbers N + , N_, the stepsize of the system fl26l) are Ax = 1/N+, Ay = 1/N_. 
So the smaller of motor numbers N + , AL, the easier for motor numbers n + , n_ to jump 
from one of the steady subdomains I, II or III to another. Intuitively, the probability 
that (n + /N + , n_/N_) lies in the neighborhood of the stable steady state point M, is 
proportional to the area of Mj's steady subdomain. Mathematically, the probability 
of motor numbers n+, n_ change from n^t , n_ to n+ , along trajectory S is 



n 

{Si,Si-\.\)&Sn 

n 

(S fc ,S fe+ i)GSj7 



n 



t 3 + 



(Sj,Sj+i)£SL 



n j+ + e j+ + TTj_ + Ci . 



7r fc+ + e fc+ + 7r fc _ + e fe . 



n 



(Si,Si + i)eSn 



(29) 
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where S L U S R U S v U S D = S, (Pi,P 2 ) G S R if and only if n+(P 2 ) = n+(Pi) + 
l,n_(P 2 ) = n_(Pi), (Pi,P 2 ) G S L if and only if n + (P 2 ) = n+(Pi) - l,n_(P 2 ) = 
n_(P t ), (P l5 P 2 ) G ^ if and only if n+(P 2 ) = n+(Pi),n_(P 2 ) = n_(P x ) + l, (P 1; P 2 ) G 
S'd if and only if n + (P 2 ) = n+(Pi), n_(P 2 ) = n_(Pi) — 1. So, theoretically, we can ob- 
tain the probability that motor numbers n+, n_ change from the neighborhood of one 
stable steady states to the neighborhood of another stable steady states. From these 
transition probabilities, we can know more details about the steady state movement 
of the cargo in this small iV + , iV_ cases. 

6 Conclusion and remarks 

In this paper, the steady state properties of the recent tug-of-war model, which 
is provided by Lipowsky et al to model the movement of cargo, which is transported 
by two motor species in the cell, is discussed. Biophysically, the stable steady states 
are the most important states, because the transition time to the stable steady state, 
as illustrated in this paper, is very short (see Figs. [7] and so almost all of the 
data are measured in stable steady states. Through the discussion in this paper, we 
can know that the final steady states of the movement of the cargo is determined by 
initial numbers of the plus and minus motors which are bounded to the microtubule. 
Certainly, the velocity and direction of the movement are also determined by other 
several parameters, such as N±, F s ±, n±, eo±, Fd±, vf±, vb±, Fext, 7- One can also find 
that, almost each of the parameters has critical points, which determine the stable 
steady velocity and direction of the cargo. It is most probable that, many of the 
parameters, including the numbers iV + and iV_ of plus and minus motors which are 
tightly attached to the cargo, and the initial binding numbers n + (0) and n_(0), can 
be determined by the biochemical environment and properties of the cargos, so some 
of which can be transported from the plus end to the minus end, and others can be 
transported reversely. 

Acknowledgments 



16 



This work was funded by National Natural Science Foundation of China (Grant 
No. 10701029). The author thanks professor Hong Qian of University of Washington 
for his help to complete this research. The author also thanks the reviewers for their 
help to improve the quality of this paper. 

References 

[1] David D. Hackney. Processive motor movement. Science, 316:58-59, 2007. 

[2] N. J. Carter and R. A. Cross. Mechanics of the kinesin step. Nature, 435:308-312, 
2005. 

[3] Y. Taniguchi, M. Nishiyama, Y. Ishhi, and T. Yanagida. Entropy rectifies the 
brownian step of kinesin. Nature Chemical Biology, 1:342-347, 2005. 

[4] R. D. Vale. The molecular motor toolbox for intracellular transport. Cell, 
112:467-480, 2003. 

[5] H. Sakakibara, H. Kojima, Y. Sakai, E. Katayama, and K. Oiwa. Inner-arm 
dynein c of chlamydomonas flagella is a single-headed processive motor. Nature, 
400:596-589, 1999. 

[6] A. M. Hooft, E. J. Maki, K. K. Cox, and J. E. Baker. An accelerated state of 
myosin-based actin motility. Biochemistry, 46:3513-3520, 2007. 

[7] J. Christof, M. Gebhardt, Anabel E.-M. Clemen, Johann Jaud, and Matthias 
Rief. Myosin- v is a mechanical ratchet. Proceedings of the National Academy of 
Sciences of the United States of America, 103:8680-8685, 2006. 

[8] Katsuyuki Shiroguchi and Kazuhiko Kinosita Jr. Myosin v walks by lever brow- 
nian motion. Science, 316:1208-1212, 2007. 

[9] H. Noji, R. Yasuda, M. Yoshida, and Jr. K. Kinosita. Direct observation of the 
rotation of fl-atpase. Nature, 386:299-302, 1997. 



17 



[10] J. Howard. Mechanics of Motor Proteins and the Cytoskeleton. Sinauer Asso- 
ciates, Sunderland, MA, 2001. 

[11] Hongyun Wang and George Oster. Energy transduction in the fl motor ofatp 
synthase. Nature, 396:279-280, 1998. 

[12] Masayoshi Nishiyama, Hideo Higuchi, and Toshio Yanagida. Chemomechanical 
coupling of the forward and backward steps of single kinesin molecules. Nature 
Cell Biology, 4:790-797, 2002. 

[13] Hongyun Wang. Chemical and mechanical efficiencies of molecular motors and 
implications for motor mechanisms. Journal of Physics: Condensed Matter, 
17:S3997-S4014, 2005. 

[14] Yunxin Zhang. The efficiency of the molecular motors. Journal of Statistical 
Physics (in press), 2009. 

[15] K. Svoboda and S.M. Block. Force and velocity measured for single kinesin 
molecules. Cell, 77:773-784, 1994. 

[16] Peter Reimann. Brownian motors: noisy transport far from equilibrium. Physics 
Reports, 361:57, 2002. 

[17] Hongyun Wang. Several issues in modeling molecular motors. Journal of Com- 
putational and Theoretical Nanoscience, 5:1-35, 2008. 

[18] Hong Qian. The mathematical theory of molecular motor movement and chemo- 
mechanical energy transduction. Journal of Mathematical Chemistry, 27:219- 
234, 2000. 

[19] P. Reimann, C. Van den Broeck, P. Hanggi H. Linke, J.M. Rubi, and A. Perez- 
Madrid. Giant acceleration of free diffusion by use of tilted periodic potentials. 
Physical Review Letters, 87:010602, 2001. 

[20] R.D. Astumian. Thermodynamics and kinetics of a brownian motor. Science, 
276:917-922, 1997. 



[21] R Dean Astumian. Biasing the random walk of a molecular motor. Journal of 
Physics: Condensed Matter, 17:S3753-S3766, 2005. 

[22] Steffen Liepelt and Reinhard Lipowsky. Kinesins network of chemomechanical 
motor cycles. Physical Review Letters, 98(25) :258102, 2007. 

[23] Michael E. Fisher and Anatoly B. Kolomeisky. Simple mechanochemistry de- 
scribes the dynamics of kinesin molecules. Proceedings of the National Academy 
of Sciences of the United States of America, 98 (14): 7748-7753, 2001. 

[24] B. Derrida, J. L. Lebowitz, and E. R. Speer. Shock profiles for the asymmetric 
simple exclusion process in one dimension. Journal of Statistical Physics, 89:135- 
167, 1997. 

[25] Steven P. Gross, Michael A. Welte, Steven M. Block, and Eric F. Wieschaus. Co- 
ordination of opposite-polarity microtubule motors. The Journal of Cell Biology, 
156:715-724, 2002. 

[26] Sean W. Deacon, Anna S. Serpinskaya, Patricia S. Vaughan, Monica Lopez Fa- 
narraga, Isabelle Vernos, Kevin T. Vaughan, and Vladimir I. Gelfand. Dynactin 
is required for bidirectional organelle transport. The Journal of Cell Biology, 
160:297-301, 2003. 

[27] Michael A. Welte, Silvia Cermelli, John Griner, Arturo Viera, Yi Guo, Dae- 
Hwan Kim, Joseph G. Gindhart, and Steven P. Gross. Regulation of lipid-droplet 
transport by the perilipin homolog lsd2. Current Biology, 15:1266-1275, 2005. 

[28] G. A. Smith, L. Pomeranz S. P. Gross, and L. W. Enquist. Local modulation of 
plus-end transport targets herpesvirus entry and egress in sensory axons. Pro- 
ceedings of the National Academy of Sciences of the United States of America, 
101(45):16034-16039, 2004. 

[29] Melanie J.I. Miiller, Stefan Klumpp, and Reinhard Lipowsky. Tug-of-war as 
a cooperative mechanism for bidirectional cargo transport by molecular motors. 



19 



Proceedings of the National Academy of Sciences of the United States of America, 
105(12):4609-4614, 2008. 

[30] Janina Beeg, Stefan Klumpp, Rumiana Dimova, Ruben Serral Gratia, Eberhard 
Unger, , and Reinhard Lipowsky. Transport of beads by several kinesin motors. 
Biophysical Journal, 94:532-541, 2008. 

[31] Stefan Klumpp and Reinhard Lipowsky. Cooperative cargo transport by several 
molecular motors. Proceedings of the National Academy of Sciences of the United 
States of America, 102:17284-17289, 2005. 

[32] Melanie J.I. Miiller, Janina Beeg, Rumiana Dimova, Stefan Klumpp, and Rein- 



hard Lipowsky. Traffic by small teams of molecular motors. arXiv:0807.0964vl 
[cond-mat.stat-mechj '. 

[33] Frank Jiilicher and Jacques Prost. Cooperative molecular motors. Physical Re- 
view Letters, 75:2618-2621, 1995. 

[34] K. Svoboda and S.M. Block. Force and velocity measured for single kinesin 
molecules. Cell, 77:773-784, 1994. 



20 



